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Abstract 

The contemporary cosmic expansion is considered in the context of Modified Friedmann Dy- 
namics (MOFD). We discuss some relativistic model exploring analogy to MOND modification 
of Newtonian dynamics. We argue that MOFD cosmologies can explain fraction of dark matter 
in the accelerating Universe. We discuss some observational constraints on possible evolutional 
MOFD scenarios of cosmological models coming from SN la distant supernovae. We show that 
Modified Newtonian Dynamics can be obtained as a Newtonian limit of more general relativistic 
models, with polytropic component of Equation of State. They constitute a special subclass of 
generalized Cardassian models basing on generalization of the Raychaudhuri equation rather than 
on generalization of the Friedmann first integral. We demonstrate that MOND cosmologies are 
compatible with observed accelerated phase of expansion of current universe only for high value of 
cosmological constant. The Bayesian framework of model selection favored this model over ACDM 
model if mjO is fixed but this evidence is not significant. Moreover obtained from statistical anal- 
ysis value of the MOND characteristic j3 parameter is far from value required for explanation of 
the flat rotation curves of spiral galaxies. 
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I. INTRODUCTION 

The idea that dark matter manifestation in flat rotation curves of spiral galaxies [3j is 
a consequence of Modified Newtonian Dynamics (MOND) seems to be both intriguing and 
controversial how claims Lue and Starkman [l[ U\- In principle there is a simple way to 
test this theory 3| by consideration observational consequences coming from cosmological 
models basing on these modifications of gravitational interactions at the late time. Of course 
the MOND theory is Newtonian but it should be obtained as a limit of more general and 
fundamental theory (Modified General Relativity - MOGR). If we consider homogeneous 
and isotropic cosmology in this theory then the modified Friedmann equations will describe 
the evolution of the Universe. 

Recently, different modifications of Friedmann equation were proposed \&\. Motivation 
of this model is to explain the current acceleration of the Universe without references 
to the unknown form of dark energy. In this scenario (called the Cardassian or polytropic 
expansion) there is no dark energy component but the universe is matter dominated. It is 
accelerating due to the adding a certain additional term to Friedmann first integral of Ein- 
stein equations (with the Robertson- Walker symmetry). An important question is whether 
there is any connections between the MOND driven cosmology and the Cardassian one. 

In the paper by Lue and Starkman [1[ it is presented interesting idea of derivation mod- 
ified theory of gravity from constraint coming from the fundamental Birkhoff law. In this 
approach authors explain how the cosmic acceleration is generated through these modifica- 
tions. 

In this paper we incorporate MOND for the late time cosmological scenario, while early 
stages of evolution are dominated by usual matter described in a standard way by general 
relativity. For simplicity (without losing the generality) we assume that our universe is flat. 
In derivations of basic dynamical equations both for Newtonian and relativistic model we 
use particle-like description 7]]. In this approach the evolution of the universe is represented 
by a motion of a unit mass particle under the action of a one-dimensional potential V(a) 
which can be simply obtained from the MOND gravitational acceleration postulate. The 
position variable a is a scale factor of the universe and all dynamics is determined by the 
potential function through the analog of Newtonian equations. 

The Cardassian models base on the generalization of the Friedmann first integral by 
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adding in r.h.s. a term which is called the Cardassian term, i.e., 

H 2 = P - + Bp- - * (1) 



where p is the energy density. If the source of gravity is a perfect fluid with pressure p = 7P 
(7 = const) then p = p (a/ao)~ 3( - 1+ "'\ a® is a present value of the scale factor, k — 0, ±1 is 
the curvature constant, H = (In a)' is the Hubble function. Note that equation (DD) is the 
first integral of the generalized Einstein equation for the Robertson- Walker symmetry. The 
basic equations constitute the system 



# = -# 2 -§(l + 3 7 ) + ^r, (2a) 
P = -3#p(l + 7) (2b) 



with the first integral in the form 



TT \ 2 / n \ - 3 ( 1+ T) / n \ -Mi+l) / n \ -2 



— = fi 7 , — + fi C ard,0 — + fifc.O ~ ■ (3) 

HoJ \a J \a J \a J 

Equation ff2a| is called the Raychaudhuri equation, while equation pbl) is the conservation 
equation. It seems to be more natural to generalize the Raychaudhuri equation instead 
of its first integral. In the last case we obtain a more general theory containing MOND 
cosmoloe; special case. If we consider the standard Cardassian models, then in the 

right-hand side of equation [3]) only power low terms can of type can appears while in the 
MOND cosmologies some part of potential is logarithmic type. 

Our basic idea is to explain the fraction of dark matter in the Universe in analogy to the 
Milgrom [2I explanation of flat rotation curves of spiral galaxies, i.e. in terms of the MOND 
conception rather than mysterious dark energy. 

The organization of our paper is as follows. In section 2 we provide a brief summary of the 
features of the Cardassian models and generalized Cardassian models that are relevant for 
our further discussion. The particle-like description of MOFD cosmologies and constraining 
model parameters in the light of SNIa data (based on the Riess sample) are presented in 
section 3. Finally in section 4 some concluding remarks and perspectives for analysis of 
cosmology in the new MOFD (or MOGR) paradigm are formulated. 
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II. GENERALIZED CARDASSIAN MODELS AS A NATURAL GENERALIZA- 
TION OF FRW MODELS. 



By the generalized FRW Cardassian models we understand models which dynamics is 
governed by the generalized Raychaudhuri equation and conservation condition 



a 2 \3 v 'J 

p = -3^(l + 7(a))p(a), (4b) 

where a dot denotes differentiation with respect to the cosmological time t, f(p(a)) defines 
the type of modification of the standard Raychaudhuri equation for FRW cosmology which 
holds for / = 0. 

System (T4J) has a first integral in the form 

a 2 k 

PeS ~ = 3^j, (5) 
er or 

where p e s(a) plays the role of effective energy density (see Appendix). Equation §5§ is 
independent (directly) on the special form of matter the Universe is filled with. In the 
generic case if we put into (j4"aj) f(p) oc p n and 7(a) like for a mixture of noninteracting 
matter and radiation then the usual class of the Cardassian models is recovered. However 
let us note that ([5]) with p c fr(a) = p(a) + 3Bp n does not play the role of the first integral 
in the special case when f(p)a oc a -1 . It is just the case of the MOND cosmologies. To 
illustrate this let us consider the simplest case of single fluid with energy density p and 
7 = const. Then from equation (14b |) we obtain 

\ -3(1+7) 

or in term of density parameter = 

/ a \ " 3(1+7) 

a = «„„(-) (7) 
Let us substitute 7(p) = 3Bp n . Hence ( l4a|) assumes the form 

It would be useful to consider in ([8]) two cases 

, 2 2 

n j± — or n = — ^. 

^3(1 + 7 3(1 + 7) 



In the first case we obtain a class of cosmologies called the Cardassian models. They can 
be treated as standard cosmological models where the universe is filled with a mixture 
of non-interacting perfect fluids with the equation of state p = 7p for the first and p = 
(71(1 + 7) — l]p = wp for the second one. Therefore the Cardassian models with a single fluid 
have dynamics in the form of a two-dimensional dynamical system 

x = y (9) 

y = -\ {fi 7l0 x- 2 - 37 (i + 3 7 ) + n CavdfiX - 2 - 3w (i + 3w)} = -dV/dx, (10) 

where f2 7>0 = p 7i o/3ifo, ficard,o = B P™,o/ Hq(1+3w) , are density parameters of matter and the 
fictitious Cardassian fluid respectively, x = a/a^ is a dimensionless scale factor in the units 
of its present value ao- A dot here denotes differentiation with respect to re-parameterized 
time variable r defined as t —> r : dt \H \ = dr. Of course fi 7j0 + f2 Cardi0 + fi fc ,o = 1 is satisfied. 
Note also that the potential function V is determined modulo to any additive constant. 
Because n ^ z(x+i) s y s t em f flO]) possesses the first integral in the form 

£ + V(x)=0, (11) 

where V(x) = — \ {fi 7i ox -1-37 + ficard,o^~ 1_3li ' + ^fc,o} and the constant in V should be 
chosen such that Yli ^«,o = 1- The constraint condition ^ ^i,o = 1 in the general relativity 
reveals the fact that both matter and geometrical term are related. The fact that ficard,o 
does not contribute in this relation is a reflection of the fact that we are beyond the standard 
cosmology. Note that it can be estimated only from the observations. 

Let us comment now the second case of n = 2/3(1 + 7). Then the potential function 
assumes very special form with logarithmic component: 

V(x) = ~ {fi 7 ,o£-^ 37 + ^c,o lnx + fi fc ,o + (1 - ^ 7 ,o)} , (12) 

where f2| ff = Qk,o + (1 — ^ 7 ,o)- Then if we substitute this form into thelfTTT) we can obtain 
the form of the first integral for this case (see Appendix). Note that both last two terms in 
( TT2|) of the same type can be defined in one term which we called effective curvature density 
parameter. Usually the form of first integral (fTTj) is treated as a starting point to further 
analysis of the generalized Friedmann equation. In our opinion the generalization of FRW 
equations on the level of the Raychaudhuri equation seems to be methodologically more 
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correct procedure than generalization of its first integral. Moreover is more general because 
one additional case is included. 

Finally the generalized Cardassian models in our terminology constitute larger class of 
models and both cases for which both n ^ 3(i+ 7 ) an d n — 3(i+ 7 ) belongs to this class. There 
are two parameters characterizing models of this class (n, fi 7 ,o) if ^fc,o = 0. 



III. PARTICLE-LIKE DESCRIPTION OF MOND AND MOFD COSMOLOGIES. 

In MOND the gravitational acceleration g exerted by a body of mass M at the radial 
distance a obeys the relationship 

J -a~ 2 for \g\ > g 
gcc* (13) 
^ -a for \g\ < g 

where go is a critical value of acceleration. Hence the potential of the gravitational field can 
be simply calculated from the formula 

m Jo I In a for \g\ < g 

jasing on the particle-like 



At first we can build the Newtonian (fTBl) cosmological models 

description of quintessential cosmology developed by us earlier [tj]. Following this approach 
the dynamics of Newtonian cosmological models can be represented by a motion of the 
particle-universe under the action of a one-dimensional potential V = V(a), a is a scale 
factor of the universe plays the role of positional variable. 

The heuristic method of obtaining Newtonian modified potential is basing on considera- 
tion Schwarzchild solution of relativistic model. We start from Newtonian model potential 
and then derive relativistic model. In the FRW cosmology the evolution of the universe can 
be derived from the Hamiltonian which in terms of dimensionless variable takes the form 

K = \y 2 + V(x) = 0, (15) 

where x = a/a ,y = x and V(x) is in the form 

V(x) = -~ {Oy.ox" 1 - 37 + ^mond.o lnx + (1 - fi 7 , )} . (16) 
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The last term (1 — f2 7) o) in (jT6j) plays only the role of negative curvature term f2£ ff = 1— £l Tj o. 
The Hamiltonian is defined on zero energy level H. — E — 0. The motion in the configuration 
space is defined in the domain admissible for motion: 

V = {x: V{x) < 0} (17) 

From ( fT5l) and (fT6l) we obtain the counterpart of the Friedmann equation in our theory. Of 
course if we substitute H = H and x — 1 then we recover Q i = 1 as a constraint on 
density parameters from relation 

H\x) = Hi {n 7i0 aT 3(1+7) + ^mond^" 2 In a; + fi c fc >~ 2 } (18) 

or in the terms of redshift 

H\z) = Hi {fi 7>0 (l + z) 3(1+7) - fi MO ND,o(l + zf In (1 + z) + Q| ff (l + z) 2 } . (19) 



By comparing (118j) with (1121) we find strictly correspondence between a special second class 
of the Cardassian models with n = „,, 2 , and MOND cosmologies. 

3(1+7) to 

From the Newtonian analogue of the equation of motion x, = —dV/dx we find that 

X = l ~ {fi 7 ,o(-l " 3 7 )X~ 2 " 37 + fiMOND.QX" 1 } (20) 

The universe is accelerating at the present epoch (x=l) if only 

^MOND.O > (1 + 37)1X^0 (21) 

Therefore for 7 = (dust) ^mond,o > ^m,o is required if f2£ ff = 0. 

The values of model parameter (f2 m o, ^mond,o) can be obtained from the fitting procedure 

to SNIa data. The luminosity distance as a function of redshift is given in the form 

, , , l + z f z dz' 
dAz) 



Ho Jo ^Q m ,o{l + z'f - fi MO ND,o(l + z') 2 In (1 + V) + (1 - f2 m , )(l + z') 2 ' 

(22) 

In the mentioned before paper by Lue and Starkman l|, |2| we find very interesting idea 
of the MOND law of gravitational interacting derived from the general relativity. The 
authors assuming the validity of the Birkhoff theorem and derive the basic cosmological 
model equation in the form 

H 2 x 2 ( p \ I + C\Qm 3 for the Einstein regime fi m > Q c 

Hi x 2 9 V Peri J ~ \ (3tt 2 J?\ntt m + C 2 fim 3 for the MOND regime Q m < fi c 
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In the Lue and Starkman 



2| model the evolution of the universe consists of two phases, 
the first one dominated by gravity following the general relativity and the second one by 
its modification (MOGR). Our idea is little different because we assume that both effects 
are acting as different regimal effects but general relativity dominates at early stages of 
evolution while MOGR describes the late time evolution. The dependence of the Hubble 
function describes following formula 

H 2 (x) = Hi {n mfi x~ 3 + ^mond.o^- 2 Ins + ^ x~ 2 } , (24) 

where fij^o is the effective curvature such that 

tfg = P£%flfotom,o + to kt0 , (25) 

Q^o = 1 - On,0, (26) 
Q«r, 

'm,0' 



^mond,o — — 3/3f2^o- (27) 



Finally we obtain the same governing equation as ([161) from the general relativistic consid- 
erations. 

The dynamics of this model can be represented in the form of the autonomous dynamical 
system 

x = y 



where 



(28) 
y = —dV/dx 



V(x) = -- {fl^ox- 1 + f] M OND,ohix + fig,} . (29) 



The system ( 1281) has the first integral in the form 

H = \y 2 + V(x) = 
From the first integral we obtain (Il{z) relation 

d L {z) = — — / - (30) 

H ° J ° Jfi m ,o(l + ^) 3 -^MOND,o(l + ^) 2 ln(l + 2 / ) + ^fo( 1 + ^) 2 



The parameter (3 can be expressed as a function of f£ m ,o- For example for dust matter 7 = 
we obtain constraint on (3 parameter 

p=- \; s Qm '° ■ (si) 



TABLE I: Results of the statistical analysis of the MOND model obtained from the best fit with 
minimum \ 2 - F denotes fixed value of parameter. 



sample 


k,0 




' *A/fni\rn n 






M 


x 2 


P 


Gold 


-0.88 


0.00 


-2.00 


1.88 


15 


.935 


173.1 


oc 






0.00 


-0.81 


1.00 


15 


.955 


175.2 


oo 




-0.90 


0.01F 


-2.00 


1.89 


15 


.935 


173.1 14.36 






0.01F 


-0.78 


0.99 


15 


.955 


175.2 


5.60 




-0.94 


0.05F 


-1.95 


1.89 


15 


.935 


173.1 


4.78 






0.05F 


-0.68 


0.95 


15 


.955 


175.3 


1.67 




-0.83 


0.30F 


-1.19 


1.53 


15 


.955 


173.8 


0.88 






0.30F 


-0.02 


0.70 


15 


.955 


175.8 


0.01 




0.00 


1.00 


1.80 




15 


.965 


177.6 


-0.60 






1.00 


1.80 




15 


.965 


177.6 


-0.60 



If we define moment in the evolution of the universe at which both material and 



MOND terms are equal we obtain 



1 + ^eq 



^MOND,o/^m.O — ; — T, — ; \j (32) 

ln(l + z eq ) 

where ft M OND,oMn,o = -3/3fi~^ 3 - 

The results of our analysis are based on the Gold Riess Riess et al. [8j) supernovae la 
sample and there are presented in the table [B On can see that considered models well fited 
SNIa data. However MOND model required value of (3 ~ 15 for possibility of explanation 
of flat rotation curve. We obtained such value of (3 only for the model with low fl m ,o — 0.01 
and J7|q = —0.9 (i.e Q^q = —3.97. This value of fi m0 and £\ i0 are in disagreement with 
both results of CMBR and primordial nucleosynthesis. 



IV. MOFD MODEL VERSUS ACDM MODEL IN THE LIGHT OF BAYESIAN 
INFORMATION CRITERION. 

In this section we extended previous model by adding dark energy in the form cosmological 
constant or phantoms We show that the MOFD cosmologies can be obtained as a Newtonian 
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limit of class of Phantom models which base on a simple modification of the FRW equation. 
The physical status of both MOFD and Phantom models is similar because they offer the 
possibility of alternative explanation of dark matter and dark energy, respectively. We 
investigate some observational constraints on the FRW cosmological models with baryonic 
matter and MOFD phase squeezed in the evolutional scenario between the epoch of matter 
domination and the dark energy epoch. We compare such a model with the concordance 
ACDM model and argue that while both models are indistinguishable (close value of \ 2 ) the 
Akaike and Bayesian informative criterions favors MOFD model with baryonic dark matter. 

We consider two possible model with the exit on A epoch or on the phantom (Cardassian) 
epoch. For both cases the relations H(z) are (respectively): 

H = #0^,0(1 + zf + Wl + z) 2 ~ fiMOND,o(l + zflnil + z) + n A (33) 

and 

H = H Qy /si mfi (l + zf + fi, cfr , (l + zf - fi MO ND,o(l + z)Hn(l + z) + n Phfi (l + z)^ (34) 

where p = wp, w < — 1 for phantoms, n = 1 + w for Cardassian. 

To compare considered models, how they fitted the data the informative criteria can be 
useful {9]. The problem of classification of the cosmological models on the light of information 
criteria on the base of the astronomical data was discussed in our previous papers jlO. 11. 

The Akaike information criterion (AIC) is defined in the following way 

AIC = -21n£ + 2rf (35) 

where £ is the maximum likelihood and d is the number of the model parameters. The best 
model with a parameter set providing the preferred fit to the data is that minimizes the 
AIC. 

The Bayesian information criterion (BIC) introduced by Schwarz is defined as 

BIC = -2\nC + d\nN (36) 

where N is the number of data points used in the fit. 

This criterion gives a simple objective criterion for the inclusion of new parameters into 
the standard ACDM model. From the results presented in the Tables HTllI VI we can draw the 
following conclusion. The Qmond is needed as a parameter and hence it is more likely that 
observations were generated in MOFD. 
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case 


name of model 




H{z) 




free parameters 


d 





Einstein-de Sitter 


H 


= #0V^m,0(l + ^) 3 




Ho, fl m ,o 


2 


1 


ACDM 


H 


= # \AWi + ^) 3 


+ n fc)0 (i + z) 2 + n A 


Ho, ^ m ,05 


3 


2a 


MOND, n m>0 - fitted 


H 




+ QkeSfii 1 + Z ) 2 ~ ^MOND,o(l + z) 2 ln(l + z) 


+ $7 a ^m,0> ^MOND,0 


n A 4 


2b 


MOND, n mfi = 0.05 








Hq, ^MOND,0, ^A 


3 



TABLE II: The Hubble function versus redshift for analyzed scenarios. 



case 


AIC - ^a,o 


= 0) AIC (l-n^o - ^a,o 


^ 0) BIC (1-O m>0 - ^a, = 0) BIC 


(l-fi^O " O A ,0 + 0) 





325.5 


194.4 


328.6 


200.5 


1 


179.9 


179.9 


186.0 


189.0 


2a 


181.2 


181.1 


190.3 


193.4 


2b 


179.3 


179.1 


185.4 


188.3 



TABLE III: The values of AIC and BIC for distinguished models (Table HU). 



V. CONCLUSIONS 



Our general conclusion is that MOND cosmology should be treated as a potential alter- 
native to the ACDM model in the context of explanation of dark matter. To clarify the 
status of these model let us consider two sets of best-fitted (gold sample of SNIa, Riess et 
al. 8J) model parameters (see tabled]). 

Our point of view is following - because the ACDM model fits SNIa data as well as 
the MOND alternative and additionally the second model explain dark matter content in 
term of ^mond,o the model under consideration should be treated as a possible candidate to 
explain dark matter in the Universe. However in particular case for the flat universe with 
^m,o — 0.3, Q\ t o = 0.7 we obtain ^mond.o — —0.02, i.e. (3 = 0.01 while for the flat universe 
with f2 m>0 = 0.05, Qa,o = 0.95 we obtain Omond^ = —0.68 i.e. (3 = 1.67. The first case 
is corresponding to the ACDM model while the second should be treated as an alternative 
description of acceleration driven by cosmological constant and dumping by baryonic matter 
(^m,o = 0.05). Both models are indistinguishable — close values of x 2 ( see table HI) and as 
result overlapping Hubble diagrams. 

On can see that however considered models well fited SNIa data MOND model required 
value of (3 ~ 15 for possibility of explanation of flat rotation curve. We obtained such value 
of j3 only for the model with f2 m0 = 0.01 and fi^o = —0.9 (i.e Qk,o = —3.97). This value 
of Q m o and Q^ are in disagreement both with result of CMBR and early nucleosynthesis. 
Finally we conclude that, the MOND conception explain only separately flat rotation curves 
of spiral galaxies or the fraction of dark matter in the Universe but it is not able to explain 
these both facts together. 

In this paper we also demonstrate that classical MOND conception can be derived from 
more fundamental relativistic theory, namely from the generalized Cardassian model. 

The main aim of the paper was to show that the existence of the MOND phase during 
the evolution of the Universe, before the epoch of domination of dark energy can explain 
the presence of dark matter in the Universe. In other words there are two indistinguishable 
scenario from the point of view of explanation of the SNIa data. On the other hand if 
we a priori assume that Q m> o — 0.3 the observations exclude the cosmological model with 
the squeezing MOND phase in the cosmological scenario (^mond,o — 0). If we assume flat 
universe with the value of f2 mj o ~ 0.3 as it is suggested by extragalactic observations than we 
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obtain that ^mond,o should be small, but not necessary equal to zero. In our approach we 
check whether the MOND phase frozen in the cosmological scenario according to Starkman's 
idea can give us understanding of the fraction of nonbarionic matter in f2 m - We find that 
such a model well fit supernovae data but value of (3 is far from Lue and Starkman Q value 
P = 15 

The second topic of this paper is the construction of the new class of cosmological models 
with frozen the MOND phase into evolutional scenario with exit to the Cardassian models. 
As it is well known the Cardassian models are an alternative to the cosmological models 
with dark energy in the explanation of present acceleration of the current Universe. In these 
models instead of dark energy violating the strong energy condition is postulated a simple 
modification of the Friedmann first integral. In this paper the model is fitted to observations 
of distant SNIa using the Riess sample. We obtain analogous results as in the case with 
the exit to the dark energy epoch. The advantage of the model with frozen MOND phase 
and exit to the Cardassian models is twofold. First, it can explain the acceleration of the 
Universe. Second, it can explain the fraction of the dark matter. 

The other results can be summarized as follows. We propose the theoretical description of 
cosmology MOFD based on the modified gravity. We find the connection of such models with 
recently discussed Cardassian models [^(J . The parameter (3 characterizing the MOND phase 
is estimated. We also estimated this parameter for the model with exit to the Cardassian 
model. In this case we obtain the value of characteristic parameter (3 which is far to the 
value assumed by Starkman (j3 = 15). The value of Cardassian exponent in the term p n 
in the modified Friedmann equation is close to zero. This situation is very close the model 
with the cosmological constant but nevertheless n is negative and nonzero. 

In this paper we pay attention to the flat cosmological models. It would useful to make 
some remarks on the non-flat cosmological models. In the models with exit to A epoch we 
estimate the curvature type term fi| ff . From the x 2 analysis we obtain that non-flat case is 
more preferable than its flat counterpart. The similar dependence of cIl(z) on the Hubble 
diagram is obtained for fitting the model with the exit to the Cardassian domination epoch 
with the SNIa data (without any prior on n) 

However MOFD cosmologies are compatible with observed late-time accelerated expan- 
sion of contemporary universe. The popular method of apriorical generalization of Fried- 
mann equation is adding polytropic component of r.h.s. of H 2 relation i.e. generalization 
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Friedmann first equation. Our proposal is generalization Raychaudhuri equation rather 
then Friedmann. Then we obtain previous generalization plus one exceptional case which is 
strictly related with main subject of the paper. 

however we still share the opinion expressed by Sahni that there is the fundamental dif- 
ficulty of MOND gravity because this theory is not embedded within a more comprehensive 
and fundamental theory of gravitation. We also do not know the Lagrangian for the Car- 
dassian modification of gravity but these models can be treated as a simple modification of 
the cosmological models with FRW symmetry. 
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VII. APPENDIX 



In this section we demonstrate how the presence of additional term in the Raychaudhuri 
equation can be modeled by some noninteracting fictitious fluid X with energy density px{a) 
and pressure px(a). We start from the basic equations 

- = -^(p + 3p) + fa m , (37a) 
a o o 

p = -3H(p + p). (37b) 

If we postulate that 

-\(px + Zpx) = ^a m (38) 



then (1371) can be rewritten to the form 

- = ~\ J2(Pk + 3p k ), (39) 
a 6 ^-^ 

i,X 

where the summation should be performed over all components of fluid. For any % fluid 
conservation equation is satisfied 

Pi = -3H(pi + pi). (40) 

Of course analogical condition should be satisfied by the fluid X, i.e 

dpx 3 



da a 
1(3 



(px+Px). (41) 



From (138]) we calculate px and then we substitute this expression into (|4TI) . Hence we obtain 

Px = -\ P x ~ ja m (42) 



and 



dpx(a) 2 B m 

-px{a)-\ a . (43) 



da a a 

As a solution of (143]) we obtain 



4 + -r^a m for m ^ -2 



Px(a) = { a m+ (44) 
-% + 4lna for m = -2 



-A - l^a™ for m ^ -2 

3a^ 3 m+2 



p x (a) = < ^ am + 2 ' (45) 

b lna f m = _ 2 

Of course system (1371) has the first integral in the form 

(X k ^ ^ . . (X , . 

i 

where 

p i = -3H(p i +p i ) (47) 

for any z-fluid, ^ pi = P an d also p G ff = —3H(p c s +p e s)- The first integral fj4"6l) has different 
form for both distinguished cases 

p + ^ + — ^-a m -3^ = 3^, for m ^ -2 (48a) 

C B a k 

p+ — + — lna - 3— = 3— , for m = -2 (48b) 
ar ar a z er 

We require the correspondence with standard FRW model for the case 5 = 0. Hence we 
obtain C = 0. Finally the potential functions for both cases takes the following form 

Up+ -^a rn )a 2 for m ^ -2 
V(a) = { eKH m+2 J r (49) 

-|(p + In a)a 2 for m = — 2 

Of course the Hamiltonian system is still determined on the energy level E = —k/2 {H = 
a 2 /2 + V(a) =E). 
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